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Abstract 

A detailed calculation is presented of the dielectric function in superconduct- 
tors consisting of two Josephson coupled superconducting layers per unit cell, 
taking into account the effect of finite compressibility of the electron fluid. 
From the model it follows, that two longitudinal, and one transverse optical 
Josephson plasma resonance exist in these materials, for electric field polar- 
ization perpendicular to the planes. The latter mode appears as a resonance 
in the transverse dielectric function, and it couples directly to the electrical 
field vector of infrared radiation. A shift of all plasma frequencies, and a re- 
duction of the intensity of the transverse optical Josephson plasmon is shown 
to result from the finite compressibility of the electron fluid. 



1 



Typeset using REVT^ 



I. INTRODUCTION 



In recent years the 'second' Josephson plasma resonance (JPR) phenomenon has been 
studied theoretically and has been observed experimentally [1-6] in double layer high Tc 
cuprates. Originally the prediction of a transverse optical Josephson resonance and the 
associated second JPR in layered superconductors with two layers per unit cell appeared as 
a short conference paper [2] , explaining briefly the main theoretical ingredients and results. 
One of the reasons for revisiting this problem is, that experimentally the transverse optical 
JPR was observed [7,8] in SmLao.8Sro.2Cu04_5 with an oscillator strength much smaller than 
expected on the basis of the simple expression derived in Ref . [2] . Following a suggestion by 
Bulaevskii [9], the expressions for the dielectric function will be re-derived in this manuscript, 
while taking into account an extra term representing the finite compressibility of the electron 
fluid. The central result is, that the dielectric function is 



which is the same expression as in Ref. [2] , except that the volume fractions xj and xk have 
been replaced with effective weight factors zj and 5/^ = 1 — 5/, which depend on the volume 
fractions, the plasma frequencies uji and ujki and the compressibility. 

Let us consider the optical response function for a material with two superconducting layers 
per unit cell. The lattice constant along the c-direction, perpendicular to the planes, is d. 
The layers are grouped in pairs, with interlayer distances xid and xxd alternating {xi-\-xk — 
1). Hence the z-coordinate of the m'th plane normalized by the lattice constant is Xm — fn/2 
if m is even, and Xm = (m — l)/2 + if m is odd. 

The charge fluctuation of each plane is characterized by a charge amplitude Qm and a phase 
0m, where m is the layer index, and d is the length of the unit cell along the c-direction. 
The discussion in this paper will be restricted to the charge fluctuations perpendicular to 
the planes, corresponding to a homogeneous charge distribution within each plane. In this 
case the electric flelds are perpendicular to the planes. The electric held of a single plane 
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with an area A, and charged with positive charge Q is constant in space, and the field fines 
are directed away from tfie plane. As a result tfie potential energy of a positive test cfiarge 
witfi cfiarge q in tfiis field decreases linearly as a function of tfie distance d from tfie plane: 
Epat{z) — —AnqQz/A. Tfie potential energy stored in tfie field of tfie cfiarge fluctuations 
Qm and Qn of a pair of planes at a distance d\xm — Xn\ is fialf tfiis amount 

Vc[Qm, Qn) = TTp:, (J-j 

wfiere Co = A/ (And) is tfie capacitance of two planes at a distance corresponding to tfie 
lattice parameter d. Tfie total potential energy stored in tfie fields of tfie cfiarge fluctuations 
is just tfie linear superposition of tfie contributions from all pairs of planes in tfie crystal. 
In tfiis context it is important to point out, tfiat tfie planes considered fiere are truly two- 
dimensional (2D) in tfie electrodynamical sense: Because tfie cfiarge fias no spatial degrees 
of freedom perpendicular to tfie planes (otfier tfian tunneling between planes) , tfie individual 
planes provide no cfiannel for metallic screening for fields along tfie c-axis. Tfiis is quite dif- 
ferent from tfie situation encountered in classical Josepfison junctions between tfiick metallic 
layers. In tfie latter case individual planes do screen tfie electric fields polarized along tfie 
c-direction. 

Tfie second source of potential energy is tfie tfie electronic compressibility. Tfiis is due to 
tfie fact tfiat if SN electrons are added to a plane, tfie free energy increases witfi an amount 
SF — jjiSN -\- SN"^ /{2Kqv?), wfiere ji, Kq, and n are tfie cfiemical potential, tfie electronic 
compressibility, and tfie electron density respectively. For a Fermi-gas K^n^ = dn/dfi cor- 
responds to density of states at tfie Fermi level. In tfie context of 'excitons' in two-band 
superconductors, tfie compressibility term fias been first considered in 1966 by Leggett [10]. 
In neutral fluids tfie compressibility causes propagation of sound, wfiereas for electrons it 
causes tfie dispersion of plasmons. Plasma dispersion of tfie JPR in tfie cuprates fias been 
described by Koyama and Tacfiiki [11]. In part tfie compressibility can be motivated by 
calculations based on tfie random phase approximation(RPA), sfiowing a finite dispersion 
of tfie cfiarge density fluctuations in single layer [12] and bilayer [13] cuprates, and of spin. 



amplitude and phase collective modes [14]. It has been shown by Artemenko and Kobel'kov, 
that the frequency of the resonance, its dispersion and its damping are strongly influenced 
by the presence of quasi-particles at finite temperature [15]. Although this type of calcula- 
tion demonstrates, that the Pauli exclusion principle causes a a finite compressibility of the 
electron fluid, leading to finite dispersion of the collective modes, weak coupling approaches 
are not well supported due to the strong electronic correlations in these materials. We the- 
fore treat the electronic compressibility as a phenomenological parameter in this paper. 
The terms in the free energy proportional to 5N only shift the equilibrium density. In har- 
monic approximation, the charge fluctuations around equilibrium follow from the quadratic 
terms 

For later convenience the dimensionless constant 70 = l/{A7rde^Kon'^), proportional to the 
2D bulkmodulus, is introduced here, to characterize the compressibility. In Ref. [2] we left 
these terms out of consideration. 

We calculate the longitudinal dielectric function following the usual procedure of adding 
external charges to each layer, distributed such as to provide an external electric field 
D (the displacement field) of the plane wave form with wavevector k. The definition of the 
dielectric constant [16] implies, that the internal and external charge distributions interact 
only via the electromagnetic field. Hence, the interaction between internal and external 
charge is described by the Coulomb term, Eq. 1, but the external charge does not enter the 
compressibility term, Eq. 2. The charge dynamics enters via the Josephson couphng J^"^^ 
between each set of nearest neighbor planes 

Hkin = -^m^^ COS(0^ - (pm+l) (3) 

Our aim is to determine the dielectric constant, and collective modes in the absence of 
external DC magnetic fields. For this purpose we will need the equations of motion for 

the internal charge accelaration d^Q\^/dt^ subject to the fields of the internal and external 
charges = + Qm- These follow from the Hamiltonian 
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m>n ^<^0 m ^<^0 m 

Here the phases (pm and the internal charges are conjugate variables, which are subject 
to the Hamilton- Josephson equations of motion: ^dcpm/dt = dH/dQl^, and ^dQ\^/dt — 
—dH/d(f)jm where e* = 2e is the charge of a Cooper-pair. 

II. EQUATIONS OF MOTION 

Working in the linear response regime, we obtain the equations of motion for the internal 
charge accelaration d'^Q]^/dt^ subject to the fields of the internal anc? external charges 

Due to the lattice periodicity, the solutions have to be plane waves with a wavevec- 
tor k = (p/d. Therefor wc can use generalized charge coordinates i?, and 5* defined as 
Q2m — 'S'e*'"'^, Q2m+i — -Rc*'"'^, to describe the charge fluctuations in the even and odd 
planes. The Josephson coupling energies J|^"^^ = / and J|^_i = K characterize the two 
types of junctions. 

It is quite easy to extend this to the situation, where we have a lattice polarizability char- 
acterized by a dielectric constant e] and for each type of Josephson junction. This 
corresponds to the transformation Co Cav, 7o 1, xj —>■ zi, and xk zk in Eqs. 
4, and 5. Here the following deflnitions have been used: Cav — ^e^^/47rci for the average 
capacitance, 7 = / {^T^de^ K^n^) for the compressibility, zi — xie^^/ej and zk — xr^Iv/^k 
(together satisfying zk + zj = 1) for the weight factors, and l/e^^ = xk/(^k xj/ej for the 
average dielectric constant. We are now ready to formulate the equations of motion for the 
generalized coordinates S and R for each wavenumber k — (f)/d. 
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2Cav {^} S' = -T.n e'^^{\n\{I + K) - \n - zi\I - \n + zk\K} 

-j:ne'^'"{\n + zi\{I + K) - \n\I - \n + 1\K} + 2'y {I + K)S' - 2'y {I + e-''^K)R' 

(6) 

2Ca. {^Y R' = - En e"^''{\n\{I + K)-\n + zj\I - \n - Zk\K} R' 

- E„ e^"^" {\n - zi\{I + K) - \n\I - \n - 1\K} + 27(7 + K)R' - 2^{I + e''^K)S' 

The convergent lattice sums over n can be replaced with the identities 

V |n + p ± xle^"-^ = e-^^'^ ^ / ^ , (7) 

n 1-COS0 

Before we continue it is convenient to partly diagonalize Eq. 6, by transforming the expres- 
sions to new generalized charge coordinates Q and P, defined as Q = S+R and P = e^'^S+R, 
with identical transformations for the internal and external charges. We will see later, that Q 
and P correspond to charge fluctuations across the barriers of type K and I respectively. The 
reverse transformations are S = {Q — P) / {1 — e^"^) , and R = [P — e^'^Q)/ {l — e^'^). In addition 
it will turn out to be convenient to introduce the Josephson plasma resonance frequencies 
characteristic of the two types of junctions — ZKK{e* /TifC~^ , and to} — zil{e* /UYC~^ , 
and the corresponding local charge response functions ex = 1 and e/ = 1 —uj'j/uj'^. 

With the help of the identities Eq. 7 the equations of motion of the external and total 
generalized charge coordinates become 

= Ci^g* - 2-iZK^uj'j^uj-^{Q' - P'e-'^'^cos (0/2)) 

(8) 

= e/P*-2727^L.'fLj-2(p«_gv'^/2cos((/>/2)) 

We see, that if the compressibility term 7 = 0, these equations of motion are already 
diagonal, corresponding to two non-dispersing longitudinal plasmons at frequencies uik and 
cu/. It is also immediately clear from this, that Q, and P correspond to the charge fluctuations 
across junctions of type and / respectively. If 7 is flnite, the equations of motion are 
coupled, and the plasma frequencies will have a flnite dispersion as a function oi k — (p/d. 
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III. CALCULATION OF THE DIELECTRIC FUNCTION 



We are interested in the response of the total electric field E to an external field D which 
is polarized along the c-axis, and which varies harmonically in time and space along the 
c-direction, i.e. 

D{f, t) = Doze'^'''-^*') (9) 
The dielectric function e{k,u!) is calculated from the definition D — eE. We will employ the 

— * — * — * — * 

fact, that D — WV^ and E — WV. V, V"^, and are the total, internal and the external 
voltages respectively. Thus we need to arrange the external charges in a such a way, as 
to garantee that the external field D{z) has a plane wave form. For z coinciding with the 
coordinates of the conducting planes, this requires that V'2m+i/^2m — ^^^''^ satisfy this 
requirement by giving the external charges the ratio R^/S^ = {e^^''^ — zk — Zie^'^}/{1 — 
zkg^^''^ — Zie~^^'^'^}. For the generalized coordinates this implies 



ZKsm{zi(j)/2) 



(10) 



ziS\n{zK(t>l'2) 

For the calculation of e we need to calculate the macroscopic average of and corre- 
sponding to the macroscopic electric and displacement fields. For this it is sufficient to know 
the voltages at the positions of the planes. Now, using Eq. 1 we observe, that the voltages 
in each plane are 

V2m ^Y.n{Q2n\n-m\+Q2n+i\n-m + Zi\) ^^^^ 

We consider a charge oscillation with wavevector k — The charges in the alternating 
layers are — e^^'^S'* and Qln+i — e"*'^-R*. After summation over the index n, the voltages 
in the zeroth and first plane 

^0 = c:j^){S^ + R\^K + zre-^^)} ^^^^ 
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with similar expressions for the external charges,for which Cav must be replaced with Co- 
The electric fields integrated between the nearest neighbor planes are 

with similar expressions for D, for which Cav must be replaced with Cq. In the limit A; ^ 
the macroscopic electric field is just the sum of the two integrals devided by the lattice 
parameter d. We conclude from this, that for A; — > the dielectric function is 



I^D{z)dz ^ zkQ' + ziP' 
!^E{z)dz ^""^ zkQ' + ziP' 



From Eq.lO we see, that P^ = Q'^ in the limit k ^ 0. We can combine this identity with 
the equations of motion, Eq. 8, to prove that P*/Q* — {ex — 2^{u;j( / zk + ujj / zi)u;~'^} / {ej — 

2^{uj],/zK+ujj/zi)uj-^}, that P7P* - €i-{{€K-ei)2^Luj/zi}/{eKOj''-2^iiul/zK+ioj/zi)}, 
and Q'^/Q^ = eK — {{^i — ^k)2ju!x/ zk}/{^ii^'^ — 2j{u!j^/zK+i^'j/zi)}. The dielectric function 
is now easily obtained, 

1 1 (cj2 _ Cj^) 



(15) 



with the definitions 



2 _ /'I I 7 ^,,,2 I /'I I 1 \, ,2 ^ . I ( (I ^ 1 \,.,2 (1 _i_ 7 V ,2^^ I i^^^K^l)'^ 



= (1 + + (I + - V((i + - (i + + (16) 

Here we have adopted the convention for labeling the two plasma-resonances such, that K 
always refers to the highest plasma resonance frequency. 

IV. CENTRAL RESULT 

Provided that ui < Cut < Cuk, it is always possible to express ujt as a weighted average 
of the two longitudinal frequencies 
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ojrj, = zkoJi + ziQj^ (17) 

with weight factors satisfying zk ^ zi = \. The latter are no longer the volume fractions zj 
and zki as in Ref. [2]. Instead they depend on the volume fractions, andora. the microscopic 
electronic parameters characterizing the two types of junctions. The effective fractions can 
be calculated by inverting the above relation, i.e. using 

ZK = {ilpK-^T)l{^\-~^''l) (18) 

As a result [17] we can write the inverse dielectric function, Eq.l5 as a linear superposition 



of two plasma resonances 



~ 2 ~ 2 

Cat) , ^K^ 



.2 ' , ,2 r,2 (19) 



which is the same expression as in Ref. [2], except that zj and zk have been replaced with 
effective volume fractions 



2 2(2:^^2;, + 27 + 472) ^ (2;j^2;, + 27)(cui-cu|)2 zkZi + 2^ + Cj\ - 

(20) 

Together, Eqs. 19 and 20 form the central result of this paper. Prom the intensities of the 
experimental loss functions the effective volume fractions zk and zi can be extracted. These 
can be used to calculate 7, and this in turn can be used to determine the density of states 

- = (21) 

or, using Ki^r? = dn/dfi, the compressibility. Once a)/ and ojk have been measured, and 7 
has been calculated from the weight factors z~k and zj using Eq. 20, it becomes possible to 
make two further deductions, namely the determination of cuk and cui using 



col = \ ± 



ZK 






I ^ 






Zj+27 1 


2 





/ 4(27)2^2^^2 








2 V 





(22) 



from which we can calculate directly the Josephson coupling energies 

2 r . -1 2 



^"d7=|-{^} (23) 
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V. EVOLUTION OF THE OSCILLATOR STRENGHT AS A FUNCTION OF 7 



For the analysis of experimental data Eqs. 19 and 20 suffice to deduce the microscopic 
parameters /, K, and Kqu^, i.e. the two Josephson energies and the compressibility factor. 
To predict the plasma resonance frequencies ujt, <^i, and ujk, we can use Eq. 16. The 
intensities of the two peaks in the energy loss function, Im— l/e(a;), are just the weight 
factors z~K and zj. Their dependence on the microscopic parameters cuj, ujk, and Kqu^ is 
given by the expressions 

1 {zk — Zi){u:]^ — lJ]) — 2^:{^'j^/ zk + / zj) , . 

Zk = i-- Zi= - ^ ; ^ {^V 

The compressibility is characterized by the dimensionles parameter 7. The most important 
effect of this extra term is, that the intensity of the highest plasma resonance, ujk, is reduced 
compared to what it would have been if 7 were zero: It is clear from Eq. 24, that the effective 
volume fraction z~k is smaller than zk- 

The oscillator strength of the transverse optical plasmon follows directly from the pole- 
strength of the pole near ut in From Eq. 15 it follows that St equals (a;|- — (jJt){^t ~ 
(^j)/(^T- With the help of this, it follows that the dependence of St on the microscopic 
parameters cuj, cuk, and Kon^ is 

{ujI - U} f {ZIZK + 27) 



(25) 



In Fig. 1 the plasmon strengths z/, and zk, and the the pole-strength St are displayed as 
a function of 7 for a few different sets of ujki oji, and the volume fraction zk- 



VI. DISPERSION OF THE PLASMA MODES 

In the previous section we found for the double layer cuprates two longitudinal plasma 
modes with the electric field vector polarized along the c-direction. The dispersion of the 
longitudinal modes can be calculated, by realizing that for longitudinal modes D = by 
virtue of the fact that longitudinally polarized free photons don't exist. Hence in Eq. 8 the 
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external charge coordinates — — 0. The corresponding 2 x 2-matrix is easily solved, 
providing the two longitudinal branches 



i - i\ + i)-'K + (I + + V (d + ^Wk - (I + + cos^ ¥ (26) 

'^f = (i + + (l + - y ((I + -{\ + i Vzf + 

An example of this dispersion is given in Fig. 3. In Fig. 4 the oscillations of these modes are 
sketched. The longitudinal modes in the lefthand side of the figure have a finite dispersion, 
provided that 7 7^ 0, in other words, if the electron gas has finite compressibility. 

In addition to the two longitudinal plasma modes there is one transverse optical plasma 
mode, and one transverse polarized pole of the dielectric function at a; = representing the 
superconducting dielectric response for fields and currents polarized along the c-axis, and 
with the direction of propagation parallel to the planes. In the righthand panel of Fig. 4 the 
oscillations of these modes are sketched. The transverse mode is coupled to electromagnetic 
radiation for long wavelengths, giving rise to coupled plasma-polariton modes instead of 
separate photons and transverse optical plasma-modes. The plasma-polariton dispersion 
follows from Maxwell's equations in dielectric media, and is given by the relation between 
wavenumber k (in the solid) and frequency (a;) 

= e(cu)a;2 (27) 
With the dielectric constant given by Eq. 15 we get 



1/2 



2 



(28) 



This dispersion is sketched in Fig. 5. We see, that there are two plasma-polariton branches. 
The lowest starts at frequency Uj in the long wavelength limit (small k), and quickly merges 
with the transverse optical plasma frequency cut, as the wavelength is reduced below c/lut, 
which is of the order of a mm if lot/^tt is of order 300 GHz. At these shorter wavelengths the 
character is almost purely the TO-JPR. The upper branch corresponds to a conventional 
transverse Josephson plasmon [18] (without the adjective 'optical') as it merges with the 
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light-line in the short wavelength limit. The lower branch is novel. It corresponds to a real 
polarization wave, and these modes can be used to convert electromagnetic radiation into 
microscopic currents, or vice versa. 

Finally it should be added, that in the short wavelength limit, the presence of finite com- 
pressibility gives rise to non-trivial dependence of e{k, u) on the wavevector k, which for the 
transverse modes is parallel to the planes. This effect becomes prominent on short wave- 
length's of the order of lattice spacings. Fig. 2 was sketched on the scale of k of the order 
of a few inverse millimeter. On the scale of a few inverse A the k-dependence of e{k, ui) will 
give rise to finite k-dispersion of u;i{k), ujT{k) , and CjK{k). 

VII. CONCLUSIONS 

An expression was derived for the dielectric function of superconducttors consisting of 
two Josephson coupled superconducting layers per unit cell, taking into account the effect 
of finite compressibility of the electron fiuid. In this model two longitudinal one transverse 
optical Josephson plasma resonance exist. The latter mode appears as a resonance in the 
transverse dielectric function, and it couples directly to the electric field vector of infrared 
radiation. A shift of all plasma frequencies, and a reduction of the intensity of the transverse 
optical Josephson plasmon is shown to result from the finite compressibility of the electron 
fiuid. 
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FIG. 1. Oscillator strength's as a function of 7. 
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FIG. 2. Oscillator strength's as a function of 7 
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FIG. 3. Dispersion of the longitudinal JPR's. 
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FIG. 4. Snapshot of the currents (arrows) and planar charge fluctuation amphtudes (indicated 
by gray-scales) of the two sets of transverse and longitudinal modes with polarization along the 
c-direction. On the righthand side of each plot the voltage distribution is indicated. 
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FIG. 5. Dispersion of the transverse JPR's for small values of fey (about one millionth of the 
Brillouinzone) . 
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